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We study the possibility of realizing quantum anomalous Hall effect (QAHE) with tunable Chern 
number through doping magnetic elements in the multi-layer topological insulator film. We find 
that high Chern number QAHE phases exist in the ideal neutral samples and it can make transition 
to another QAHE phase directly by means of tuning the exchange field strength or sample thickness. 
With the help of an extended Haldane model, we demonstrate the physical mechanism of the tunable 
Chern number QAHE phase. Further, we study the phenomena in the realistic samples with the 
Fermi energy being gating-modulated, where the multiple-integer QAHE plateaus remain and two 
separated QAHE phase are connected by metallic phase. The stability of the high Chern number 
QAHE phase is also discussed. 

PACS numbers: 73.43.Nq, 72.15.Rn, 72.25.-b, 85.75.-d 



I. INTRODUCTION 



The integer quantum Hall effect!^] is one of the 
most important discoveries in condensed matter physics. 
When a strong perpendicular magnetic field is applied 
to the two-dimensional electron gas at low temperatures, 
the Hall conductance exhibits a precise quantization in 
units of e 2 /h (i.e. the fundamental conductance unit) 
due to Landau-level quantization. This quantization is 
directly connected to the topological properties of the 
two-dimensional bulk states, characterized by a topolog- 
ical invariant C known as the first Chern number! 3 ^ 

Non-zero Chern number can in principle also occur in 
the band structure of other systems with time-reversal 
symmetry breaking, such as in a ferromagnet, leading 
to the so-called quantum anomalous Hall effect (QAHE). 
This effect was first proposed in a seminal paper by Hal- 
dane^ in a honeycomb lattice model with an average 
flux per unit cell being zero. Due to its unique non- 
trivial topological properties and great potential appli- 
cation for designing dissipationlcss spintronics, extensive 
studies have been made to search possible host materi- 
als of realizing such QAHE. Several candidate systems 
are proposed recently including the ferromagnetic Mer- 
cury based quantum wells in the insulating stated the 
disorder induced Anderson insulator P Rash ba gr aphene 
coupled with exchange fieldP Kagome latticepE2l anc j op- 
tical lattice modelsP-S 

However, although the QAHE has been theoretically 
proposed, the Chern numbers in the proposed QAHE 
systems can only take some low and limited values (i.e. 
C = 1, 2, 4). This is distinct from that in the con- 
ventional integer quantum Hall effect, where the Chern 
number can be very high and tuned continuously through 
adjusting the external magnetic field or the Fermi en- 
ergy. Therefore, a natural question arises: whether it 
is possible to find some host materials exhibiting QAHE 




FIG. 1: (Color online) Schematic plot of the setup: a mag- 
netically doped topological insulator on a substrate with its 
Fermi energy being tuned by the attached external gate. Blue 
arrows indicate the magnetized dopants. N z measures the 
sample layers along Z-direction. 

with tunable Chern numbers? In the following, we shall 
give a positive answer by doping magnetic elements into 
the topological insulator films. Topological insulator is 
new quantum phase of matter that behaves as an in- 
sulator in the bulk but carries time-reversal symmetry 
protected odd-number pairs of helical edge states (two- 
dimensional) or surface states (three-dimensional) 

The topological insulator was first proposed in two- 
dimensional materials, however, there are only limited 
candidates showing the theoretical possibility of hosting 
two-dimensional topological insulator i n contrast, 
there have been numerous materials or compounds found 
to host the three-dimensional topological insulators both 
theoretically EIHSD anc j experimentall}! 2 ^!!. This has at- 
tracted extensive attentions from the condensed-matter- 
physics community, and provides an interesting plat- 
form for the practical application in low-power dissipa- 
tion spintronics. For instance, it has reported that mag- 
netically doped topological insulator could exhibit the 
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QAHE using first principal theory,^ and recent exper- 
imental discovery has shown the giant anomalous Hall 
conductance in the magnetic doped topological insulator, 
which is a good sign for the final realization of QAHE. 31 

In this article, we show that the QAHE with tunable 
Chern numbers can be achieved through doping mag- 
netic elements in a topological insulator film. The mag- 
netic dopan t in the topological insulator results in two 
effectsP^ESI ( a ) breaking the time-reversal symmetry; (b) 
inducing an effective exchange field. Through investigat- 
ing the evolution of the bulk band structure, counting the 
edge states winding number and the direct Hall conduc- 
tance calculation, we determine the existence of QAHE 
phases with tunable Chern numbers in our studied sys- 
tem. Using an extended Haldane model, we give an ex- 
planation on the mechanism of the tunable Chern num- 
ber QAHE phases. Furthermore, the influence of Fermi 
energy fixed by external gate are discussed. When the 
Fermi energy is slightly deviated from the zero energy, 
the high Chern number QAHE phases can still survive 
even though the metallic phases emerge between two sep- 
arated QAHE phases. In the end, we study the disorder 
effect on the obtained QAHE phase and find that the 
high Chern number QAHE plateaus can be very stable 
against weak disorders. 

The remaining of this paper is organized as follows. 
In Sec. [TTJ we introduce an theoretical model and the 
methods for calculating kinds of topological features. 
In Sec. |III| we present the numerical results for the 
ideal neutral samples and discuss their mechanisms. In 
Sec. |IV| we show the influence of the Fermi energy shift- 
ing and disorder on tunable Chern number QAHE which 
are inevitably present in realistic samples. Finally, a brief 
discussion and conclusion are given in Sec. [V] 

II. MODEL AND METHODS 



as: 

a a 
T a = ^(Jo <g> T z + ^CT ® T - % A^ a a ® T x 

mo = ma z ® r (2) 

where M, A a , B a , D a are independent parameters with 
M determining the inverted band gap amplitude and A a 
reflecting the Fermi velocity, a is the unit vector along 
a = (x, y, z) direction, a and r are the Pauli matrices in 
spin and orbital spaces, respectively. Hi mp is used to de- 
scribe the magnetic dopants, m measures the effective ex- 
change field strength. For simplicity, we assume that the 
Lande g factor is the same for all bases, which is reliable 
because in realistic three-dimensional topological insula- 
tor all lowest-energy bands determining the band topol- 
ogy are always combined by the p-orbitals^or d-orbitals. 
Since we focus on the resulting phenomena of the topo- 
logical insulator film, the layers along the Z-direction are 
set to be finite. 

In general, analyzing the evolution of the bulk and 
edge energy spectrum as functions of some tunable pa- 
rameters is an efficient way to investigate the topological 
features of a system. Let us first focus on the bulk en- 
ergy spectrum. Since the studied system has the trans- 
lational symmetry along both x and y directions, both 
the corresponding momenta k x and k y are good quan- 
tum numbers. Through performing the partial Fourier 
transformation 

= 7^^^"^^'' (3) 



Figure [T] plots the schematic setup: the topological in- 
sulator material is doped with magnetic elements, with 
the arrows representing the magnetic dopants. We use 
Nz to label the thickness of the topological insulator lay- 
ers. As a starting point, we introduce the effective Hamil- 
tonian. In the cubic lattice model, the electronic state at 
each site i can be expressed as cj>\ = [ai-f, &if, ay,, &iJ T , 
where (a, b) denote two independent orbits and (f, 4-) re P~ 
resents spin indices. In the tight-binding representation, 
the topological i nsulator doped with magnetic elements 
can be written as j 22 l 34 l 3 S l 

H = Hsd + H imp 

H3D=J2 E ^t& + J2 E 4>tTcc<t>i+a + h.C. 
i i a—x^y,z 

H unp = m 0<l>t<fii (!) 

i 

where H 3 £> describes the bulk Hamiltonian of the three- 
dimensional topological insulator, Eq, T ai mo are written 



the real-space Hamiltonian in Eq. ([T]) becomes 

ff x (k) = ^2 {<t>i xky (z)[E + m }<fi kxkv (z)} 

+ E {ti x kSz)[T*e ik * +T v e ik *)<t> kxky {z) 
+ ^ik y ( z ) T ^ kv (z + l) + h.c.}. (4) 

By directly diagonalizing the Hamiltonian in Eq. ^ , one 
can obtain the bulk band spectrum. 

In order to study the edge state physics, one has to con- 
sider a boundary or surface, i.e. in our case, we choose to 
terminate the topological insulator along y direction and 
consequently an edge exists in the y — plane. There- 
fore, only k x is left to be a good quantum number. The 
corresponding partial Fourier transformation and the re- 
sulting momentum-space Hamiltonian H2{k x ) can be ex- 
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pressed as: 

Mv,z) = Ve lfc ^0i(x,?;,z), 

v x.y 

H 2 {k x )= {<t>t {y^)[E Q +m }^(y,z)}+ ]T 

kx ,y,z k x ,y,z 

{0+ (y, z)T x e lk ^ ko: (y, z) + 0+ (y, z)T 2 fcx (y, z + 1) 
+/i.c.} + ^ {0i.(y,«)T,0 fc .(y + l,z) + /i.c.}. 

= Y^H 1D {k x ) (5) 

In this way, the system can be treated as a semi-infinite 
quasi-one dimensional tight-binding chain along y direc- 
tion. Using the nonequilibrium Green's function tech- 
nique, the edge Green's function G r (k x ,E,y — 0) can 
be numerically obtained,^ where E denotes the Fermi- 
energy. Both the bulk and edge energy spectra informa- 
tion can be included in the spectral function A(fc x , E) 
of the edge, which reads 

A(k x , E) = --lm{TY[G r (k x ,E, y = 0)]}. (6) 

In addition, the topological invariant is another impor- 
tant quantity to characterize the topological properties of 
the system. For example, in the absence of external mag- 
netic field a better way to judge QAHE is to see whether 
the Chern number is non-zero or not. The Chern number 
equals to the zero temperature Hall conductance and can 
be calculated via the Kubo-formulaP 

e 2 fdk x dk y „ (*lf£l")Hf&IQ (7 , 

where Ep is the Fermi energy, £;/„ and \l/n) are the 
corresponding eigenenergy and eigenstate of H\ , respec- 
tively. 

For simplicity, in the following numerical calculations, 
we assume that the topological insulator film is spatially 
isotropic. The independent parameters are set to be: 
A a = A = 1.5, B a = B = 1.0, D a = D = 0.1, M = 0.3. 
a denotes x, y, z and nearest hopping parameter B is 
set as the energy unit. Our main results of this paper 
still hold for anisotropic parameters or more realistic pa- 
rameters.^ 



III. QAHE IN THE IDEAL NEUTRAL SAMPLES 

In this section, we focus on the topological phenomena 
in the ideal neutral samples. Physically, in our model, 
ideal neutral means that in the clean topological insulator 
film all valence bands are occupied while all conduction 
bands are unoccupied. Moreover, magnetic doping does 
not bring in some extra carriers, and the system is not 
affected by the external environment. 




x_> r -> m X h> r M x ^> r -> M 



FIG. 2: (Color online) The bulk band spectra along the high- 
symmetry lines of the magnetic element doped topological in- 
sulator film in the presence of different exchange field strength 
m=0(a), 0.037(b), 0.3(c), 0.525(d), 0.7(e), 0.819(f), 0.9(g), 
1.089(h), 1.15(i). The layer-thickness is set to be N z = 12. 



According to the definition of the topological order, 
any two insulating states are topologically equivalent 
only when they can adiabatically change into each other. 
Therefore, the classification of insulating state is highly 
related to the evolution of bulk energy spectrum. If 
bulk gap closes and reopens, the initial state and finial 
state cannot be adiabatically connected. Consequently, 
they belong to different topological phases and undergo 
a quantum phase transition. In our system, the time- 
reversal symmetry is broken due to the presence of ex- 
change field, which is different from magnetic field in 
the conventional quantum Hall effect. In realistic ex- 
periments, the exchange field strength could be tuned 
through controlling the doping concentration of magnetic 
elements. 

In Fig. [2j we plot the bulk energy spectrum along the 
high symmetry lines for different exchange field strength 
Tji. Without the exchange field, the bulk band gap opens, 
which indicates an insulating state [see in panel (a)]. 
When the exchange field is introduced, the bulk gap grad- 
ually decreases along with the increasing of the exchange 
field to. As shown in panel (b), when the exchange field 
reaches the critical value m c = 0.037, the bulk band gap 
is completely closed. When the exchange field is further 
larger than to c = 0.037, the bulk band gap reopens [see 
in panel (c)]. This gap closing and reopening indicate a 
quantum phase transition, which has been pointed out in 
Ref. [50]. 

The major finding of this article is that the bulk 
band gap is not monotonous but oscillates as a func- 
tion of the exchange field strength to. For example, 
besides to c = 0.037, we found that there are other 
critical points (band gap closing and reopening), e.g. 
to c = 0.525, 0.819, 1.089. This indicates that there are 
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FIG. 3: (Color online) Edge spectral functions A(k x , uj) of the 
semi-infinite sample along y-direction for different exchange 
field strength m=0.0(a), 0.3(b), 0.7(c), 0.9(d). The sample 
thickness is set to be N z = 12. 
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FIG. 4: (Color online) Edge spectral functions A(k x ,uj) of 
the semi-infinite sample along y-direction for different sample 
thickness N z =6(a), 12(b), 18(c), 24(d). The exchange field 
strength is set to be m = 0.7. 



more than two topologically different phases. From the 
critical points of view, one can see that in Fig. [2] there 
are five different phases. 

One of the most striking characteristics of the non- 
trivial insulating states is the existence of topologically 
protected gapless edge states, i.e. in quantum Hall ef- 
fect the number of Chiral gapless edge states equals to 
the non-zero Chern number. Therefore, edge state anal- 
ysis can be regarded as an efficient method to reveal 
the topology of the bulk states.^ Figure [3] exhibits the 
edge spectral function A(k Xl E) for the first four topo- 
logical phases with different exchange field strength be- 
ing m=0.0(a), 0.3(b), 0.7(c), 0.9(d). Note that in all the 
figures of the edge spectral function A(k x ,E) the blue 
regime denotes bulk gap, green regime denotes bulk band 
states, and the lines represent edge states at the bound- 
ary. 

In the absence of exchange field, one can notice that 
there is no edge state inside the bulk band gap in panel 

(a) . This means that the system has the same topol- 
ogy as that of the vacuum and is considered as a trivial 
insulator. When the system enters the second topolog- 
ical phase, i.e. 0.037 < m < 0.525, one chiral gapless 
edge state appears inside the bulk band gap [see in panel 

(b) for m = 0.3], indicating that the system becomes a 
nontrivial insulator with winding number N = 1. Ac- 
cording the relationship between winding number of the 
edge state and the bulk Chern number Cp-^we claim that 
the system belongs to the C = 1 QAHE phase. After the 
system undergoes the second quantum phase transition, 
we observe that there are two chiral gapless edge states 
locating inside the bulk band gap [see in panel (c) for 
in = 0.7], which labels that our system enters another 
QAHE phase with Chern number being C — 2. By con- 
tinuously increasing the exchange field strength m, one 
can obtain N — 2,3,4,5... chiral edge states. In other 



words, the QAHE phase with various non-zero Chern 
numbers can be achieved in our system. Besides control- 
ling the exchange field amplitude, we further find that the 
consecutive topological phase transitions can also happen 
through varying the sample thickness Nz, when the ex- 
change field strength m is larger than the gap controlling 
parameter M . This can be concluded in Fig. [4] show- 
ing the evolution of the edge states with different sample 
thickness Nz. 

To give a better understanding of this QAHE phases 
with various Chern numbers, we studied their topological 
properties by directly calculating the zero-temperature 
Hall conductance. In Fig. [5j we calculated the Hall con- 
ductance a xy as functions of the exchange field m and 
the sample thickness Nz- In panel (a), the phase dia- 
gram of the Hall conductance a xy in (Nz, m) plane is 
plotted, and different colors are used to specify kinds 
of Hall plateaus. One can find that the Hall conduc- 
tance varies as functions of both the exchange field and 
sample thickness. To be more specific, in panel (b) we 
show the Hall conductance as function of the exchange 
field at fixed sample thickness Nz =6,9, 12; in panel 
(c) we plot the Hall conductance as a function of the 
sample thickness at fixed exchange field m=0.3, 0.5, 0.7, 
and 0.9. Since the topological invariant Chern number 
in the QAHE is identical to the Hall conductance (in 
units of e 2 /h), therefore each separated QAHE phase and 
the resulting phase boundaries can be easily determined. 
Figure [5] is the central result of this article. The ob- 
tained QAHE phases manifest the following features: (i) 
Multiple discrete QAHE phases exist in our system with 
phase transitions directly from one to another by vary- 
ing the exchange field strength or sample thickness; (ii) 
In contrast to the reported QAHE proposals, the Chern 
number in our proposal can be very high and the highest 
Chern number is comparable to the total layers Nz; (hi) 
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the Hall conductance a xy is not monotonous as a function 
of the exchange field strength m. For instance, panel (b) 
can be divided into three regions where Hall conductance 
<j xy is stepped by 1, -2, and 1 along with the increasing 
of exchange field strength to; (iv) We emphasize that the 
relationship between a xy and to in our model resembles 
that between a xy and the external magnetic field in a 
conventional quantum Hall effect. 

So far, we have showed the existence of the tunable 
Chern number QAHE phases and demonstrated their 
topological properties in the neutral samples. In the fol- 
lowing, we will move to the physical mechanism lead- 
ing to such phenomena. C ompared to the previous two- 
dimensional QAHE models pE 8 * 30 | 3S | , our studied model is 
quasi-three dimensional. However, to satisfy the con- 
dition of requiring the vanishing wave-function at the 
regions z < and z > N z , the wave vector k z in Z- 
direction should take real discrete values with their mag- 
nitude approach k z = ^p(n = 1, 2, 3, N z ) or one imag- 
inary value 39 SI The latter case leading to the C = 1 
QAHE phase reported in Ref. J5D]- We will focus on the 
former one. In the momentum space, the Hamiltonian 
of our system can be written as H(k) = J^ fc H kz ^v) 
with k x , k y being the momenta along x and y direc- 
tions and k z taking some concrete real constants. In 
this way, our model can be regarded as the combina- 
tion of a series of two-dimensional square lattice mod- 
els described by H kz (k x ,k y ). Moreover, H kz {k x ,k y ) can 
be block-diagonalized and the new diagonal Hamiltonian 
H kz (k x ,k y ) is expressed as: 



H kz {k x ,k y ) 

h kz (k x , k y ) 
h-1 z (k x , k y ) 



£k z {k Xl k y )+^ k * v ^ ^2 



O 

{k x , ky) 



d\a z 
dl<j z 



d x & X 



d y a y 
d y a y 



(8) 



where the parameters d xi cL, d 2 z are written as: 



^k z (fax i fcy ) 

d x 
d„ 



—3D + D(cos k x + cos k y + cos k z ) 
A sin k x 
A sin kit 



dl(k x ,k y ) = -yj [M kz (k x ,k y )] 2 + A 2 sin 2 k z + ■ 



d z {k xi k y 



[M kz (k x ,k y )} 2 + A 2 sin 2 k z 



m 



M kz (k x , k y ) = M — 3B + B(cos k x + cos k y + cos k z ) 

(9) 

One can find that both two block Hamiltonian h\ and 
h 2 are the Haldane's model that is similar to the model 
described in Ref. [35]. Only when d\{dl) can change its 
sign in the first Brillouin zone, the two bands of h\ (h 2 ) 
become inverted to result in QAHE with C = +1(— 1) 
Chern numberP^For positive exchange field strength to, 
h 2 .^ is a trivial insulator since d 2 cannot change its sign. 
In contrast, d\ can change its sign and consequently h\ 
describes the QAHE for certain positive exchange field 




FIG. 5: (Color online) (a) Phase diagram of the Hall conduc- 
tance a xy in the plane (m, N z ). (b) and (c) plot the Hall 
conductance for fixed N z =6/9/12 and m=0. 3/0.5/0.7/0. 9. 
The phenomena beyond the parameter space in panel (a) are 
also shown. 



strength to. Combining the results for h\ and h 2 , the 
total Chern number of H kz in the different zone of the 
exchange field strength m can be expressed asP^l 

< to < TO c i(fc z ) or to > m c ^{k z ) 
C = I 1 m cl (k z ) < to < m c2 (k z ) (10) 
-1 m c2 (k z ) < to < m c3 (k z ) 

with the constants to c1 , to c2 , to c3 being: 



m c i( 



m c2 (k z 



(k z ) = J{M -B + Bcosk z ) 2 +A 2 sin 2 k z 



(M -3B + B cos k z 



A 2 ; 



m c3 (k z ) = J (M -5B + B cos k z ) 2 + A 2 sin 2 k z (ll) 



The Eqs. [8 |l0| are the main mechanism of the tun- 
able Chern number QAHE phases. The obtained results 
in the numerical simulations can be explained by such 
mechanism. In the following, we will demonstrate how it 
is applied to the continuously increasing exchange field 
strength to. Since k z takes a series of discrete values, the 
constants m c i, vn c2l m c3 also take a series of discrete 
values. When to exceeds m c i, the corresponding H kz 
transitions from the trivial insulator to the C — 1 QAHE 
phase. And consequently the quantized Hall conductance 
a xy is increased to e 2 /h. a xy will keep as a constant un- 
til to exceeds another m c \. Therefore, the quantum Hall 
plateaus appear in Fig. [5j As shown in panel (b), when 
Hall conductance a xy reaches its maxima, it will begin 
to drop 2e 2 /h due to QAHE transition from C = 1 to 
C = — 1 for one H kz . For even larger to exceeding m c3 , 
H kz transitions from C = — 1 QAHE phase to a trivial 



6 



3 

_ 2 



(a) 







0.0 








-E F = 


0.05 






0.10 







(b) 



r 




0.4 0.8 1.2 0.0 0.4 0.8 




FIG. 6: (Color online) (a)-(b): The Hall conductance a xv 
as a function the exchange field strength m under different 
Fermi energies Ep — 0.00, 0.05, 0.10 at fixed sample thick- 
ness N z = 9(a) and N z = 12(b), respectively, (c)-(d): The 
Hall conductance cr xy as a function the layer thickness m un- 
der different Fermi energies Ef = 0.00, 0.05, 0.10 at fixed 
exchange field m = 0.5(a) and m = 0.7(b), respectively. 



insulator. Consequently, a xy shows quantized plateaus 
with interval of e 2 /h at the right-hand-side of panel (b). 

Next, let us explain the phenomena in the panel (c). 
According to Eqs. (10) and (11), m c \ blow the fixed m 



will contribute one Chern number. Though k z and m c \ 
vary for different sample thickness N z , the numbers of dis- 
crete m c i below m will not change until N z exceeds the 
critical thickness. This explains why the Hall conduc- 
tance also shows quantized plateaus in panel (c). From 
Eq. (11) one can obtain that the minimum of m c i ap- 



proaches the gap controlling parameter M. Therefore, in 
order to observe the higher Chern number QAH state, the 
exchange field strength should be larger than M. This 
can be verified in our numerical simulations, e.g. in panel 
(c) only the C = 1 QAHE is obtained when m is identical 
to gap parameter M. 



IV. ROBUSTNESS OF THE TUNABLE QAHE 
IN THE REALISTIC SAMPLES 



In the above section, we study the topological phe- 
nomena in the ideal neutral samples. However, in the 
realistic materials such ideal conditions can hardly be 
satisfied. In general, there are two kinds of inevitable 
impurity effects. First, the impurities may bring in some 
extra carriers that make the Fermi energy shift from the 
bulk band gap. Second, the imperfect especially nonuni- 
form doping may induce the disorder effects leading to 
the magnetic or nonmagnetic scattering. In the follow- 
ing, the influence of these two effects on the high Chern 
number QAHE will be addressed. 



We first concentrate on the case that impurities pro- 
vide some extra carriers. Fortunately, one can overcome 
this difficulty by utilizing an external gate to fix the Fermi 
energy into band gap [see as figure 1]. In experiment, 
the Fermi energy can be tuned by the gate voltage. In 
Fig. [6j the Hall conductance a xy as the function of ex- 
change field strength m or sample thickness N z for differ- 
ent fixed Fermi energies Ep is plotted. The presence of 
the quantized plateaus prove that the high Chern number 
QAHE phase can still hold to some degree. Nevertheless, 
a regime with continuously varying Hall conductance ex- 
ists between two consecutive quantized Hall plateaus, 
which indicates that the metallic phases emerge between 
two separated QAHE phases. For the fixed Fermi energy 
shift, the metallic phases will easily show up in the high 
Chern number QAHE regimes with strong exchange field 
strength m or large sample thickness N z . 

As discussed in the above section, our model can be 
considered as a series of two-dimensional Haldane model 
[see in Eq. (Jsj)] that is characterized by the wave vector 
k z . Obviously, in order to obtain the quantized Hall con- 
ductance, the Fermi energy should locate inside the bulk 
band gaps of all the Haldane models, i.e. from Eq. (8) 
and (9) all possible discrete k z should satisfy the follow- 
ing relationship: 

e kz (0,0) - |4(0,0)| <E F < e kz (0,0) + |4(0,0)| (12) 



If Eq. ( 12 ) is not satisfied, the metallic phase will appear. 



A more detailed analysis of Eq. ( 12 ) shows that in order 
to find the high Chern number QAHE in the realistic 
materials, the smaller m or thinner N z as well as the 
lower Fermi energy shift are required. This analytical 
conclusion agrees well with our numerical simulations. 

Secondly, let us discuss the disorder effect. Since dis- 
order breaks the translational symmetry, it is difficult to 
study its effect on the topological properties using energy 
band theory. However, through analyzing the resulting 
transport properties, disorder effect on the bulk topolog- 
ical properties could be estimated. In the following, we 
apply our tight-binding model to a two-terminal device 
and study its transport properties using the nonequilib- 
rium Green's function method. Here we will not intro- 
duce the calculating method in detail which could be 
found in Ref. [41]. Our considered model is illustrated 
in the inset of Fig. [7] The disorders are only considered 
in the central scattering region. The source and drain 
are perfect semi-infinite leads which will not bring any 
redundant scattering at the interfaces between the leads 
and central region. In general, there are two kinds of 
disorders existing in a real system: the white noise and 
the nonuniform doping. In our model, the white noise is 
modeled by the on-site disorder energy Wi(Jo®To with W% 
uniformly distributed in the interval of [— where 
W measures the strength of disorder. The nonuniform 
magnetic doping is modeled by the uniform background 
which produces on-site zeeman energies ma z <£> tq and 
the exchange field at different sites fluctuate leading to 
mriiCTz <g)To on each site. Here rij is uniformly distributed 
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exchange field strength m 

FIG. 7: (Color online) The longitudinal conductance osd as 
a function of the exchange field strength m under different 
white noise disorder strength W and the nonuniform doping 
factor nw combinations. The error bar is used to denote the 
conductance fluctuation 5a sd- The Fermi energy is set to 
be Ef = 0.05. Inset: the schematic plot of the two-terminal 
device. Disorders are only considered in the central scattering 
regime, which is modeled by a N x x N v x N z — 30 x 40 x 9 
cubic lattice models. 

in the range [— ^f 1 , with nw being the nonuniform 
doping factor. 

Figure [7] plots the longitudinal conductance <jsd and 
the corresponding conductance fluctuation 8a sd versus 
the exchange field strength m for different combina- 
tions (W, n w )=(0, 0), (1, 0.4), (2, 0.8), and (3, 1.2). 
In absence of disorder, a$D shows perfect quantized 
conductance plateaus (see the squared-symbol curve). 
When disorder is present, we found that around the 
regimes close to the quantum phase transition (e.g. m G 
[0.94,1.06]) the conductance is easy to be destroyed in- 
dicating that the disorder results in the metallic phase 
permitting the state being backscattered. The obtained 
metallic region is in excellent consistence with the re- 
sults from the Hall conductance calculation with the 
Fermi-energy shift. While in the central regime of each 
plateau (e.g. m £ [0.74,0.86]) the longitudinal con- 
ductance as d remains quantized with vanishing conduc- 
tance fluctuations 5a sd even when the white noise dis- 
order strength approaches W = 3 and nonuniform fac- 
tor reaches nw = 1.2. Since the existence of chiral edge 
states is the only possible mechanism to interpret the ab- 
sence of backscattering, one can conclude that the high 
Chern number QAHE phases are robust against strong 
white noise and nonuniform magnetic doping P^l 

V. DISCUSSION AND CONCLUSION 

One challenge of obtaining the high Chern number 
QAHE phase is that the exchange field strength m should 



be larger than the bulk band gap amplitude of topolog- 
ical insulator. This condition is hardly to be realized in 
Bi 2 Se 3 and many other obtained materials. However, 
with the rapid development of this fielcP^, it is highly 
possible to find some narrow-gap topological insulators 
in the near future. In addition, we want to emphasize 
that the sign of gap controlling parameter M is not im- 
portant. In other words, one can obtain the tunable high 
Chern number QAHE in semi-metal or narrow-gap nor- 
mal semiconductor though the QAHE phase arisen from 
the surface state is difficult to be realized. 

Moreover, we want to discuss the relationship between 
the QAHE and the extended Haldane model. In the main 
text, we have claimed that the extended Haldane model 
can be used to explain the emergence of the quantized 
Hall plateaus. However, in the realistic materials there 
are several features that might destroy certain charac- 
teristics of the system, thus deviates the scope of the 
extended Haldane model. For example, (i) when the 
topological insulator is doped with some magnetic impu- 
rities, Lande g factors of all the occupied orbitals might 
be slightly different, thus the mass-term in the Hamilto- 
nian has diverse values and the Hamiltonian cannot be 
deduced to the Haldane model; (ii) impurities may de- 
stroy the inversion symmetry, which leads to a much more 
complicated situation that cannot be captured by the 
Haldane model. Fortunately, in our numerical simula- 
tions we have considered this two important features, and 
our results indicate that the tunable high Chern number 
QAHE plateaus are very robust against these two kinds of 
violations. This robustness can be attributed to the fact 
that the discrete finite k z protects the system from small 
disturbance and the stability of QAHE under disorders 
should be very significant in the experimental realization. 

In summary, using the tight-binding method, we 
demonstrate that tunable Chern number QAHE can be 
in principle observed through doping magnetic elements 
in the topological insulator film. The perfect quantum 
plateaus in the (m, N z ) plane exist in the ideal neutral 
samples. The topological properties of the tunable high 
Chern number QAHE are discussed via the evolution of 
bulk and edge states as well as the direct Hall conduc- 
tance calculation. Moreover, the physical origins of the 
tunable Chern number QAHE phase are explained us- 
ing the extended Haldane model. Further, we also study 
the robustness of the tunable Chern number QAHE at 
fixed Fermi-energy and nonmagnetic/magnetic disorders. 
Our theoretical prediction should shed some light on the 
searching of high Chern number QAHE materials. 
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